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Abstract 



In the context of statistical physics, Chandrasekharan and Wiese recently introduced the fermio- 
nant Fernij., a determinant-Hke quantity where each permutation tz is weighted by —k raised to the 
number of cycles in n. We show that computing Fermj. is #P-hard under Turing reductions for any 
constant k> 2, and is ®P-hard for k = 2, even for the adjacency matrices of planar graphs. As a 
consequence, unless the polynomial hierarchy collapses, it is impossible to compute the immanant 
Imm;^ A as a function of the Young diagram A in polynomial time, even if the width of A is restricted 
to be at most 2. In particular, if Ferm2 is in P, or if Imm;^ is in P for all A of width 2, then NP c Rp 
and there are randomized polynomial-time algorithms for NP-complete problems. 

1 Introduction 

Let A be an n X n matrix. The fermionant of A, with parameter k, is defined as 



where S„ denotes the symmetric group, i.e., the group of permutations of n objects, and c(7r) denotes 
the number of cycles in n. Chandrasekharan and Wiese [10] showed that the partition functions of 
certain fermionic models in statistical physics can be written as fermionants with k = 2. This raises the 
interesting question of whether the fermionant can be computed in polynomial time. 

Since (— 1)"+'^'^") is also the parity of n, the fermionant for fc = 1 is simply the determinant which 
can of course be computed in polynomial time. But this appears to be the only value of k for which this 
is possible: 

Theorem 1. For any constant k > 2, computing Fernij, for the adjacency matrix of a planar graph is 
#P-hard under Turing reductions. 

Indeed, recent results of Goldberg and Jerrum [1141] imply that for k > 2 the fermionant is hard to 
compute even approximately. Unless NP c RP, for k> 2 there can be no "fully polynomial randomized 
approximation scheme" (FPRAS) that computes Ferm/^ to within a multiplicative error 1 + e in time 
polynomial in n and 1/e with probability at least, say, 3 /4. 



n 



FermfcA=(-l)" J](-fc)^(")]~[A, 



■i,n{i) ' 



JieS,, [=1 



For the physically interesting case fc = 2, we have a slightly weaker result. Recall that ©P is the class 
of problems of the form "is \S\ odd," where membership in S can be tested in polynomial time. 

Theorem 2. Computing Ferm2 for the adjacency matrix of a planar graph is ©P-hard. 

By Valiant and Vazirani [12311 and more generally Toda's theorem 112011 , the polynomial hierarchy II18II 
reduces to ©P under randomized pol3Tiomial-time reductions. Therefore, if Ferm2 is in P then the 
polynomial hierarchy collapses, and in particular NP c RP. 

Theorems [l] and [2] imply new hardness results for the immanant. Given an irreducible character 
of Sn, the immanant Imm;^ of a matrix A is 

n 
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In particular, taking A to be the parity or the trivial character yields the determinant and the permanent 
respectively. Since the determinant is in P but the permanent is #P-hard | |21| | , it makes sense to ask how 
the complexity of the immanant varies as A's Young diagram ranges from a single row of length n (the 
trivial representation) to a single column (the parity). Strengthening earlier results of Hartmann [15], 
Biirgisser [9] showed that the immanant is #P-hard if A is a hook or rectangle of polynomial width, 

-^1 = w,'^2 = ■■■ = ■^n-w+i = 1 or = • • • = A„/„, = w , 

where w = r2(n^) for some 5 > 0. Recently Brylinski and Brylinski [l7[| improved these results by 
showing that the immanant is #P-hard whenever two successive rows have a pol3aiomial "overhang," 
i.e., if there is an i such that A; — Ai+i = n(^n^) for some 5 > 0. The case where A has large width but 
small overhang, such as a "ziggurat" where A; = w — t + 1 and n = w(w + l)/2, remains open. 

At the other extreme, Barvinok [H and Biirgisser [8] showed that the immanant is in P if A 
is extremely close to the parity, in the sense that the leftmost column contains all but 0(1) boxes. 
Specifically, [8] gives an algorithm that computes Imm;^A in time 0(s^d^n^ log n), where 5;^ and 
are the number of standard and semistandard tableaux of shape A respectively. If the height of A is 
n — c where c = 0(1), then and are bounded above by (") -/cT = 0(n'^) and (") ("^^~^) = 0(n^'^) 
respectively, giving a polynomial-time algorithm. 

Any function of the cycle structure of a permutation is a class junction, i.e., it is invariant under 
conjugation. Since any class function is a linear combination of irreducible characters, the fermionant is 
a linear combination of immanants. If fc is a positive integer, it has nonzero contributions from Young 
diagrams whose width is k or less: 

Fenn^,A = ^d'^^hmm^TA, (1) 

A 

where A ranges over all Young diagrams with depth at most k, and d^ denotes the number of semistan- 
dard tableaux of shape A and content in {l,...,k}. To see this, first note that the class function fc'^'-"-' is 
the trace of ti's action on (C'')®" by permutation, i.e., 

7r(vi (8) V2 (8) • • • V J = v„(i) (8) v„(2) <S)---<8) v„(„) . 
By Schur duality (see e.g. [13]) the multiplicity of A in (C*")®" is d^^-*, so 

A 
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We then transform fc'^*^"^ to (— 1)"(— fc)'^*^"^ by tensoring each A with the sign representation, flipping it to 
its transpose A^, which has width at most k. 

Since there are OCn*^"^) = poly(n) Young diagrams of width k or less, for any constant k there is a 
Turing reduction from the fermionant Ferm^ to the problem of computing the immanant Imm;)^ where 
A is given as part of the input, and where A has width at most k. Then Theorems [l] and |2] give us the 
following corollary. 

Corollary 3. For any constant integer k, the problem of computing the immanant Imm;^^A as a function 
of A and A is #P-hard under Turing reductions if ?c > 3, and ®P-hard if ?c = 2, even if A is restricted to 
Young diagrams with width k or less. 

In particular, unless the polynomial hierarchy collapses, it cannot be the case that Imm;;^ is in P for ^ 
of width 2. This is somewhat surprising, since these diagrams are "close to the parity" in some sense. 

We conjecture that the fermionant is #P-hard (as opposed to just ©P-hard) when k = 2. Moreover, 
we conjecture that the immanant is #P-hard for any family of Young diagrams of depth n — n^, or 
equivalently any family with a polynomial number of boxes to the right of the first column: 

Conjecture 4. Let A(n) be any family of Young diagrams of depth n — n^ for some constant 5 > 0. Then 
Imm^^^ is #P-hard. 

Roughly speaking, this would imply that the results of []3l |8]] showing that Imm;i^ is in P if A has depth 
n - 0(1) are tight. 

2 #P-hardness from circuit partitions and the Tutte polynomial 

Proof of Theorem [I] We simply recast known facts about the circuit partition polynomial and the Tutte 
pol5aiomial. A circuit partition of G is a partition of G's edges into circuits. Let denote the number 
of circuit partitions containing t circuits; for instance, is the number of Eulerian circuits. The circuit 
partition polynomial j(G;z) is the generating function 

oo 

j(G;z) = 2r,z^ (2) 
t=i 

This polynomial was first studied by Martin II17II , with a shghtly different parametrization; see also BT] 
|l[5l[IIl[l6l[2l[25]l. 

To our knowledge, the fact that j(G; z) is #P-hard for most values of z first appeared in [12J. We 
review the proof here. Recall that the Tutte polynomial of an undirected graph G = (V, E) can be written 
as a sum over all spanning subgraphs of G, i.e., all subsets S of E, 

r(G; X, y) = ^(x - (y - i)c(SH|si-lv| _ 

SCE 

Here c(G) denotes the number of connected components in G. Similarly, c(S) denotes the number of 
connected components in the spanning subgraph (V, S), including isolated vertices. When x = y, 

T{G;x,x) = Y,ix - iy(s)+m-ciG)^ 

SCE 
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Figure 1 : Left, a planar graph G (white vertices and black edges) and its medial graph (gray vertices 
and directed edges). Right, a subset of the edges of G, and the corresponding circuit partition of G^. 

where £(S) = c(S) + |S| — |V| is the total excess of the components of S, i.e., the number of edges that 
would have to be removed to make each one a tree. 

If G is planar, then we can define its directed medial graph G^. Each vertex of G^ corresponds to an 
edge of G, edges of G^ correspond to shared vertices in G, and we orient the edges of G^ so that they 
go counterclockwise around the faces of G. Each vertex of G^ has in-degree and out-degree 2, so G^ is 
Eulerian. The following identity is due to Martin [.17] ; see also [24] , or [j2j for a review. 

;(G,;z) = z^^G^nG;z + l,z + l). (5) 

We illustrate this in Figure [l| There is a one-to-one correspondence between subsets S Q E and circuit 
partitions of G^. Let v be a vertex of G^ corresponding to an edge e of G. If e e S, the circuit partition 
"bounces off e," connecting each of v's incoming edges to the outgoing edge on the same side of e. If 
e ^ S, the partition crosses over e in each direction. The number of circuits is then c(S) + £(S), in which 
case @ yields (js). 

Vertigan [|26ll proved that the Tutte polynomial for planar graphs is #P-hard under Turing reductions, 
except on the hyperbolas (x — l)(y — 1) e {1, 2} or when (x,y) e {(1, 1), (—1, —1), (co, co*), (co*, co)} 
where co = e^"'/''. It follows that computing j(G;z) is #P-hard unless z e {0, — 2, ±l,±-/2}, and in 
particular if z < —2. 

To reduce the circuit partition polynomial to the fermionant, we simply have to turn Eulerian-style 
circuit partitions, which cover each edge once, into Hamiltonian-style ones, which cover each vertex 
once. Given a directed graph G, let Gg be its edge graph. Each vertex of G^ corresponds to an edge (u, v) 
of G, and there is an edge between (u, v) and (u\ v') if v = so that (u, v) and (u', v') could form a 
path. Then each circuit partition of G corresponds to a permutation n of the vertices of G^, and 

Ferm^Ag = j(G;-?c), 

where is the adjacency matrix of Gg. In particular, if G is planar, G^ is its medial graph, G^ g is the 
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edge graph of (which is also planar) and ^ is its adjacency matrix, then 

FemifcA^^, = (-fc)^^^) r(G; 1 - fc, 1 - fc), 
and this is #P-hard for any k>2. □ 

As alluded to above, Goldberg and Jerrum II14II showed that T{G;x,y) is inapproximable for planar 
graphs and many values of x and y, and in particular in the region x,y < —1. This implies that, unless 
NP c RP, there can be no FPRAS for Fermj, for k> 2, preventing us from computing the fermionant to 
within multiplicative error e in time polynomial in n and 1/s. 



3 ©P-hardness from Hamiltonian circuits 

Proof of Theorem^ We will prove Theorem |2] by showing that the fermionant Ferm2 can be used to 
compute the parity of the number of Hamiltonian circuits in an undirected graph of size n > 4. 

Let A be the adjacency matrix of an undirected graph G with n vertices. Each permutation 7i e S„ 
that gives a non-zero contribution to Ferm2(A) corresponds to a Hamiltonian-style circuit partition of G; 
that is, a vertex-disjoint union of undirected cycles that includes every vertex. A 2-cycle in n corresponds 
to a circuit that travels back and forth along a single edge. Any circuit of length greater than 2 can 
be oriented in two ways, so a circuit partition with C2 2-cycles and c' longer cycles corresponds to 2'^ 
permutations. 

Taken together, these permutations contribute (—1)" {—2y^{—4y to Ferm2A. This is a multiple of 8 
unless n = 2 and C2 = 1, or n = 4 and C2 = 2, or if C2 = and c' = 1. This last case corresponds to a 
Hamiltonian cycle. Thus 

1 

-Ferm2A=2#H. (6) 

Valiant ['221 showed, using a parsimonious reduction from 3-SAT to Hamiltonian Cycle, that the 
problem ©Hamiltonian Circuits of computing #H mod 2 is ©P-complete. Indeed, he showed this even 
in the case where G is planar and every vertex has degree 2 or 3. Since ([6]) gives a reduction from 
©Hamiltonian Circuits to Ferm2, this shows that Ferm2 is ©P-hard as well. □ 

More generally, Ferm^i^ is at least as hard as computing the number of Hamiltonian circuits mod k. 
But for fc > 3, our #P-hardness result is stronger. 



4 Conclusion 

Is the fermionant #P-hard for k = 2? And are immanants #P-hard even for some Young diagrams of 
width 2, as in Conjecture |4]? The reduction of Theorem [l] fails in this case, since for x = y = —1 the 
Tutte polynomial is easy. In particular, if G = (V, £) then 

FermaA^ e = (-2)'^^^ r(G; -1, -1) = (-2)^^^^ (-1)'^' (-2)'^*'"'' . 

Here dimS is the dimension of G's bicycle space — the set of functions from E to Z2 that can be written 
both as linear combinations of cycles and as linear combinations of stars — which we can determine in 
polynomial time using linear algebra. 
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However, there is no reason to think that another reduction couldn't work. After all, the number 
of Eulerian circuits of G^, is the number of spanning trees of G, which is also in P — but Brightwell and 
Winkler showed that counting Eulerian circuits is #P-hard in general [6J. 

Finally, it is interesting that our proof of #P-hardness is an indirect reduction from, say, the chromatic 
polynomial, passing through the Tutte polynomial and its ability to perform polynomial interpolation by 
decorating the graph. Thus another open question is whether there is a more direct reduction to the 
fermionant from the permanent or, say, #Hamiltonian Circuits. 
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